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This  paper  investigates  the  observability  of  relative  satellite  dynamics  when  only  relative 
angle  measurements  are  available.  It  has  been  shown  previously  that,  for  non-maneuvering 
motion,  the  Cartesian  states  are  unobservable  with  angles-only  measurements  and, 
furthermore,  none  of  the  Cartesian  states  can  be  determined  from  these  measurements.  In 
this  paper  it  will  be  shown  that  if  the  relative  state  is  characterized  instead  using  a  geometric 
parameter  set  known  as  relative  orbit  elements,  many  of  these  parameters  can  be  calculated 
directly  from  the  angle  measurements.  This  method  allows  one  to  quickly  get  an  idea  of  the 
geometry  of  a  relative  trajectory  based  only  on  angle  measurements. 


I.  Introduction 

The  primary  objective  of  a  satellite  navigation  system  is  to  accurately  determine  a  satellite’s  inertial  position  and 
velocity  in  space  through  onboard  measurements  pertaining  to  the  satellite’s  position  with  respect  to  either  the 
Earth  or  navigation  beacons.  Similarly,  a  relative  navigation  system  allows  the  determination  of  the  relative 
position  and  velocity  of  one  satellite  (a  “deputy”)  with  respect  to  a  reference  satellite  (a  “chief’).  This  type  of 
system  uses  measurements  pertaining  to  the  deputy’s  position  with  respect  to  the  chief.  The  design  of  autonomous 
relative  navigation  systems  has  been  a  popular  research  area  the  past  few  years.  Such  systems  allow  a  spacecraft 
that  is  following  a  desired  relative  trajectory  with  respect  to  the  chief  to  detect  deviations  from  the  desired  trajectory, 
so  it  can  then  compute  the  required  corrective  action. 

The  most  common  measurement  types  for  relative  navigation  systems,  as  with  inertial  navigation  systems,  are 
range  and  angles;  in  this  case,  the  range  is  the  distance  from  the  deputy  to  the  chief,  and  the  angles  are  the  relative 
azimuth  and  elevation  indicating  the  direction  from  deputy  to  chief,  with  respect  to  a  nominal  plane  (usually  the 
chiefs  orbit  plane).  However,  in  some  cases  the  deputy  may  not  possess  a  ranging  sensor,  but  can  only  measure 
relative  angles.  These  angles  (azimuth  and  elevation)  can  be  immediately  converted  into  a  line-of-sight  (LOS)  unit 
vector  from  the  deputy  to  the  chief  in  the  deputy’s  body-fixed  frame.  Then,  using  the  satellite’s  attitude 
determination  system  and  a  series  of  coordinate  rotations,  this  LOS  vector  can  be  expressed  in  whatever  coordinate 
frame  is  desired. 

It  has  been  shown  previously1  that,  for  a  relative  navigation  filter  based  on  the  Hill ’s-Clohessy- Wiltshire  (HCW) 
equations2"3,  the  Cartesian  states  characterizing  the  deputy’s  relative  motion  (i.e.,  its  relative  position  and  velocity 
components  with  respect  to  the  chief)  are  unobservable  with  angles-only  measurements,  if  neither  satellite  is 
maneuvering;  namely,  none  of  the  Cartesian  states  can  be  recovered  from  these  measurements.  However,  Ref.  1 
indicated  that  the  angle  measurements  do  allow  one  to  determine  what  “family”  the  relative  trajectory  belongs  to. 
That  is,  it  can  be  determined  whether  the  trajectory  is  stationary  with  respect  to  the  chief  or  drifting,  as  well  as 
information  about  the  shape  of  the  trajectory.  Thus,  there  ought  to  be  some  amount  of  geometric  information 
recoverable  from  angle  measurements.  In  this  paper,  the  authors  choose  to  characterize  the  relative  state  using  a 
geometric  parameter  set  known  as  relative  orbit  elements.4  These  are  based  on  the  solution  to  the  HCW  equations 
and  give  a  geometric  rather  than  a  Cartesian  representation  of  the  relative  state.  It  will  be  shown  that,  for  a  given 
trajectory,  the  values  of  several  of  the  relative  orbit  elements  (ROEs)  can  be  calculated  directly  from  the  line-of- 
sight  vectors.  Because  the  angle  measurements  are  imperfect,  the  calculated  values  will  only  be  considered  an 


1  Graduate  Research  Assistant. 

2  Research  Aerospace  Engineer,  Senior  Member  AIAA. 


1 

American  Institute  of  Aeronautics  and  Astronautics 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

AUG  2008 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2008  to  00-00-2008 

4.  TITLE  AND  SUBTITLE 

Estimating  Geometric  Aspects  of  Relative  Satellite  Motion  Using 
Angles-Only  Measurements 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Utah  State  University, Dept  of  Mechanical  &  Aerospace 

Engineering, Logan, UT, 84322 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

AIAA  Guidance,  Navigation  and  Control  Conference,  Honolulu,  Hawaii,  August  2008. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

- 1 - 

- 1 - 

16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 
ABSTRACT 

Same  as 
Report  (SAR) 


18.  NUMBER 
OF  PAGES 

17 


19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


estimate.  Although  not  a  full-blown  relative  navigation  scheme,  this  method  allows  one  to  quickly  get  an  idea  of  the 
geometry  of  a  relative  trajectory  based  only  on  angle  measurements. 


II.  Relative  Satellite  Dynamics 

Three  major  assumptions  inherent  in  the  HCW  equations  are  that  the  only  force  modeled  is  that  of  a  point  mass 
gravitational  field;  the  chief  is  in  a  circular  orbit;  and  the  distance  between  the  satellites  is  small  compared  to  their 
orbital  radius.  These  assumptions  yield  the  following  linear  time-invariant  differential  equations: 

x  -  2n  y  -  3n2  x  =  0 

y  +  2nx  =  0  (1) 

z  +  n2z  =  0 


The  equations  are  written  in  the  local-vertical,  local-horizontal  (LVLH)  coordinate  frame,  whose  origin  is  at  the 
chief  satellite.  In  these  equations,  x  is  the  component  of  the  deputy’s  position  vector  relative  to  the  chief  in  the 
radial  direction  positive  away  from  the  Earth;  y  is  the  along-track  component  positive  along  the  velocity  vector  of 
the  chief;  z  is  the  cross-track  component  perpendicular  to  the  orbital  plane  of  the  chief;  and  n  is  the  mean  motion  of 
the  chief. 

The  solution  to  Equations  (1)  is: 


x  =  -^-sin(wf)  -  (3x0  +  ^^-)cos (nt)  +  (4x0  + 
n  n  n 

2x  4  y  2x 

y  =  — —  cos  (nt)  +  (6x0  +  — -)  sin(«t)  -  (6n  x0  +  3  y0  )t - -  +  y0 

n  n  n 

zn 

z  =  —  sin(^)  +  z0  cos  (nt) 
n 

x  =  x0  cos  (nt)  +  (3  nx0  +  2  y0  )  sin  (nt) 
y  =  -2x0  sin(n/)  +  (6 nx0  +  4 y0  )  cos(nt)  -  ( 6nx0  +  3 y0 ) 
z  =  z0  cos  (nt)  -  nz0  sin(«?) 


where  xo,  yo,  etc,  are  conditions  at  some  epoch  time  to,  and  t  is  the  time  since  to .  Consider  the  following  change  of 
coordinates  from  x,y,z,x,y,z  :12 


a„  =  2,||  —  |  +  |  3x  +  2  — 

n 


yd=y-  2: 


Z 


max 


xd  =  4x  + 2  — 
n 

t 3  =  atan2(x,3nx  +  2  y) 
y  =  atari  (/?z,z)-  atan2(i,3/?x  +  2v) 


(3) 


where  ae,  xj,  yd,  ft,  zmax,  and  y  are  the  ROEs.  The  inverse  of  this  transformation  is 
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-  a 


cos  ft  +xd 


x  = 


y  =  ae  sin  /3  +  yd 

^  sinfr  +  p) 


It  was  shown  in  Ref.  2  how  the  ROEs  evolve  with  time: 


a=  a 


eO 


3  3 

yd  =  yd  o  --"Xdo*  =  yd0  -~nxd 


a. 


x  =  —  n sin  6 

2 


y  =  aen  cos  /3  ~  —  nxd 

z  =  wcos  (y  +  P) 


Xd  ~  Xd0 

P  =  Po+nt 

y  =  y0+nt 


(4) 


(5) 


These  equations  are  analogous  to  Eqns  (2)  for  X,  y,  z,  X,  y,  z  in  that  they  express  the  ROE  values  at  any  given  time 
as  a  function  of  their  initial  (epoch)  values  and  the  time  since  epoch. 


The  parameterization  of  Eqns  (4)  reveals  that  the  relative  motion  in  the  x-y  plane  of  the  deputy  with  respect  to  the 
chief  is  a  superposition  of  periodic  motion  in  x  and  y,  with  period  equal  to  that  of  the  chief  ’s  orbit,  and  secular 

3 

motion  in  y.  Essentially,  this  is  an  elliptical  path  that  is  drifting  in  the  y-direction  at  a  rate  of  -  —  TlXd  .  The 


instantaneous  center  of  the  ellipse  is  (x</,  yj).  It  has  a  semimajor  axis  of  length  ae  in  the  along-track  direction  and 
semiminor  axis  of  length  aJ2  in  the  radial  direction.  /3  is  a  parametric  angle  (i.e.  phase  angle)  indicating  the  location 
of  the  deputy  satellite  in  its  trajectory,  with  (3  =  0  corresponding  to  the  perigee  location  (the  “bottom”  of  the  ellipse). 
The  relative  motion,  if  the  elliptical  path  were  “frozen”  at  a  point  in  time,  is  depicted  in  Figure  1.  Although  the 
ellipse  is  actually  drifting,  it  has  been  frozen  in  order  to  conveniently  label  the  ROEs.  The  z  motion,  according  to 
the  HCW  model,  is  purely  sinusoidal  and  independent  of  x  and  y. 


Earth 


orbital 

velocity 

direction 


Figure  1.  Planar  Projection  of  Relative  Motion  Trajectory  with  Relative  Orbit  Elements  Labeled. 
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III.  Relationship  Between  LOS  Vectors  and  ROEs 


For  this  exercise  it  will  be  assumed  that  Line  of  Sight  (LOS)  unit  vector  has  already  been  obtained  and  has  the 
components  ix,  iy,  and  iz,  pointing  along  the  LVLH  x,  y  and  z  axes  respectively.  The  LOS  vector  can  be  expressed  as 

LOS  =  -  (6) 

r 

where  r  is  the  relative  position  vector  between  the  two  satellites.  Since 

/V  /V  ^ 

r  =  xi  +  yj  +  zk  (7) 

we  can  rewrite  the  LOS  vector  in  terms  of  the  Cartesian  relative  states: 


LOS  = 


xi  +  yj  +  zk 


V 


2  2  2 
x  +  y  +  z 


Substituting  in  Eqns  (4)  will  yield  an  expression  for  the  LOS  vector  in  terms  of  ROEs: 


(8) 


^xO 

iyO 

izo 


[xd  -  Tfcos(A))]  i  +  [yao  +  aesin(po)]j  +  [zmaxsin(il; 0)]  k _ 

-  jcos2(p0))  +  ae  (2 ydosin(Po)  -  xdcos(p0 ))  +  x\  +  yj0  +  z^sin2^ 0) 


(9) 


^xl 

iyl  ~ 
i  zl 

_ [xd  -  gf  cos(fti)]  i  +  [ydi  +  OeSin{P i)]  j  +  [zmaxsin(r{f i)]  k _ 

y/a%  (!  -  fc°s2(/?i))  +  ae  (2t/dlsm(/3i)  -  Xdcos^))  Tx%  +  (10) 

Thus,  if  we  acquire  the  LOS  vector  (i.e.  angle  measurements)  at  several  times  along  a  trajectory,  eventually  we  will 
have  enough  equations  to  solve  for  the  unknowns  (ROEs). 


IV.  Results 

Equations  (9)  and  (10)  represent  six  relationships  and  six  unknowns.  Even  though  these  equations  are  extremely 
nonlinear,  and  some  of  the  ROE’s  are  not  observable,  we  are  postulating  that  some  of  the  ROE  values  can  be  found 
either  explicitly  or  numerically. 

A.  Families  of  Solutions  given  two  LOS  Vectors 

Having  a  known  p0  and  \j/0,  greatly  simplifies  the  problem  of  solving  for  the  other  LOS  values,  p  and  \j /  are 
periodic  values  that  could  be  obtained  in  a  number  of  ways.  One  simple  way  is  to  derive  P  and  \j /  from  the  azimuth 
and  elevation  angle  histories.  P  =  0  when  the  elevation  angle  is  at  a  local  minimum.  \j/  =  0  when  the  azimuth  angle 
ascends  through  zero.  An  example  using  the  initial  conditions  below  may  be  found  in  Figure  2. 


4 

American  Institute  of  Aeronautics  and  Astronautics 


ae  =  .03  km 
xd  =  -01  km 

yd0  =  .6  km 
A,  =  07r/180  radians 

Zmax  0.01  km 

ipQ  =  On/ 180  radians 


Figure  2.  Example  of  solving  for  rp  and  |3  from  Azimuth  and  Elevation  histories. 

Equations  (9)  and  (10),  relating  the  LOS  measurements  to  the  ROE’s,  may  be  reformulated  in  terms  of  the  range 
between  the  chief  and  deputy. 


lr0 


[xd  -  tfcosiPo)]  i  +  [ ydo  +  aesin(Po )]  j  +  [zmaxsin^ 0)]  k 


Rc 


2  _  [xd~  ^cos(Ai)]  i  +  [ydi  +  aesin(j3i )]  j  +  [zmaxsin(ipi)]  k 

lfl  —  p 

Rl  (11-12) 


where 


Ro  =  \  ja2e  (  1  -  ^cos2(/?o)  )  +  ae  (2ydosin(f3Q)  -  xdcos(f30 ))  +  x2d  +  yj0  +  sin2  (ip, 0) 


(13) 
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and  similarly  for  As  a  result,  the  measurements  ir0  and  iri  correspond  to  a  single  solution  associated  with  an 
assumed  value  for  R0.  The  result  is  a  whole  family  of  solutions  dependent  on  the  value  of  R0. 

Once  an  assumed  value  for  R0  is  selected,  zmax  may  be  solved  by  noting  that 

~  Roij& 

max  sin(ipo)  (14) 

Rearranging  Equations  (11)  and  (12)  results  in  the  equation  below: 


2 fad  Roixo)  (Ro^yO  Ud)  2 (Xd  Rl^xl)  (Rl^yl  Vdl) 

cos((3 q)  sin(p0)  cos(fii)  sin((3 1)  (15) 


Solving  for  this  family  of  solutions  given  ir0,  iri,  At,  |30,  and  rp0  is  necessarily  a  numerical  process.  Note  that  only  the 
first  three  equations  are  necessary  to  constrain  the  problem,  and  the  solution  is  represented  by  the  intersection  of 
three  planes  as  seen  in  Figure  3.  As  the  assumed  value  for  R0  is  increase  the  magnitude  of  xd,  yd,  and  ae  will  increase. 


Figure  3.  Representation  of  the  solution  for  xd,  yd,  and  ae  given  an  assumed  initial  range,  two  FOS  vectors  and  the 

time  they  were  apart. 

An  example  solution  set  with  unknown  xd,  yd,  ae,  and  zmax  and  the  actual  values  below  is  shown  in  Figure  4. 

v  =  0.001085410283934  rad/sec 


„..  .  At  =  198 sec 

ae  =  .1)1  km 

Xd  =  -002  km  ~ 

'  0.028383618444956  ' 

IrO  — 

0.999358093716378 

Vdo  =  -03  km 

_  0.021857967141048  _ 

flo  =  807r/ 180  radians 

'  0.055713307466721  ' 

max  0.01  km  = 

0.995603147336654 

•0Q  =  57r/180  radians 

0.075302060957631 
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x  10 


Figure  4.  Representation  of  some  of  the  possible  trajectories  given  two  LOS  vectors,  At,  oo,  |30,  and  rp0. 


B.  Solutions  Given  Additional  Info 

In  the  case  where  ae  is  known,  it  is  a  simple  matter  to  numerically  solve  Equation  (15)  for  xd,  yd0  ,  and  zn 
long  as  ae  is  nonzero.  Numerical  solver  ambiguity  about  the  sign  of  zmax  can  be  eliminated  by  also  solving 


as 


Roizo 


R\izi 


(16) 


Similarly,  in  the  case  where  xd  is  known  then  it  is  a  simple  matter  to  numerically  solve  Equations  (15)  and  (16)  for 
ae,  yao  ,  and  zmax  as  long  as  xd  is  nonzero.  Note  that  the  parts  of  Equation  (15)  that  contain  xd  must  be  used  in  order  to 
obtain  the  correct  solution.  In  the  case  where  yd  is  known,  it  is  a  simple  matter  to  numerically  solve  Equations  (15) 
and  (16)  for  ae,  xd,  and  zmax  as  long  as  yd  is  nonzero.  Note  that  the  parts  of  Equation  (15)  that  contain  yd  must  be  used 
in  order  to  obtain  the  correct  solution.  In  the  case  where  zmax  is  known,  it  is  a  simple  matter  to  numerically  solve 
Equations  (15)  and  (16)  for  ae,  xd  ,  and  yd  as  long  as  zmax  is  nonzero. 


C.  Solutions  for  Special  Cases 

In  the  case  where  ae  =  0  than  the  following  will  be  true: 


lro 


Xd 

UdO 

tc  sinty  o) 


1 

Rq 


1r\  — 


Xd 

yd0  -  %uxdA t 


1 

R~i 
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(17) 


Ri  =  +  +4i«2Wi) 

Ro  —  Arbitrary  range  value 


If  xd=0  than  the  following  holds  true: 


1 

-fc(/3 ,) 

ir0  — 

VdO  +  <7es(/?o) 

Ro 

%r\  — 

VdO  +  aes{Pl) 

Zmax  sifl  ("00  ) 

ZmaxSin^lfli ) 

Ri  =  y«e  +ae2ydls({31)  +  y2dl  +  zliaxs2(^1) 

Ro  =  Arbitrary  range  value 


In  a  practical  sense,  this  means  that  ix  and  iy  and  iz  will  follow  the  same  track  every  orbit.  Two  examples  are  seen  in 
Figure  5.  This  will  hold  true  even  when  (3  and  \\)  are  unknown.  Note  that  the  y  axis  is  a  normalized  value  with 
respect  to  the  magnitude  of  the  unit  vector. 


ix  as  a  function  of  time  \  as  a  function  of  time 


Figure  5.  LOS  measurement  histories  when  xd  =  0. 

Unless  xd  =0,  yd  will  equal  zero  for  only  an  instant.  This  is  not  very  useful,  so  only  the  situation  when  both  xd  and 
yd=0  will  be  examined.  In  this  case,  the  following  will  hold  true: 


-fc(Po) 

1 

-fc(A) 

iro  — 

aes(p0) 

~R~lrl  = 

ZmaxSin^x/j q) 

ZmaxS'i'R'i^  1 ) 

Ri 


+  Z?naxs2('lPl) 


Ro  =  Arbitrary  range  value 
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As  in  the  case  where  xd  =0,  ix  and  iy  and  iz  will  follow  the  same  track  every  orbit.  However,  the  iy  component  will 
oscillate  about  0  as  seen  in  Figure  6.  This  will  hold  true  even  when  (3  and  rf  are  unknown. 


i  as  a  function  of  time 

X 


If  zmax  =0  then  the  following  holds: 


i 

H 

a- 

1 

»°|c? 

g 

_ 1 

1 

xd  -  <fc{P i) 

iro  — 

Vdo  +  aes(P o) 

0 

Rq 

irl  — 

Vdl  +  aes(Pl) 

0 

R\  =  y  < a2  ^!0|C2/3i^  +  ae(2ydls(pi)  -  xdc(/3i))  +x2d  +  y2dl 
Rq  =  Arbitrary  range  value 


Thus,  if  two  LOS  measurements,  not  taken  a  half  period  apart,  both  have  an  iz  component  of  0,  then  zmax  must  equal 
zero.  This  will  hold  true  even  when  \j /  and  p  are  unknown. 

D.  Analytical  Solutions  with  a  Single  LOS  Measurement 

Squaring  Eqn  (9)  and  then  multiplying  both  sides  by  the  denominator  (here,  cos(P)  is  replaced  by  c(P),  and  sin(P) 
by  s(p)),  results  in: 
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1  -  -/m 


+  ae  (2 yds((3)  -  xdc((3))  +  x2d  +  y2d  +  z^s2 {$)  = 


*max*  KWJ 


(21) 


This  is  a  form  that  allows  analytical  solutions  for  many  of  the  ROE’s. 

In  the  special  case  where  a e  is  the  only  unknown  and  all  measurements  are  perfect,  it  is  possible  to  analytically 
solve  for  ae  given  any  component  of  a  single  LOS  vector  (ir  =  [ix  iy  iz  ]T  ),  but  using  all  three  often  eliminates 
multiple  positive  and  real  solutions.  This  is  done  as  follows.  Collecting  like  terms  results  in  the  following  quadratic 
equations: 


I**  (l  -  ~  °e  +  {**(2j/ds(/?))  +  xdc((3)}ae  +  ... 

{i2x(Xd  +  Vd+  Zmaxs2W)  -  Xd}  =  0 

(X  ~~  -  g2(^)}  °e  +  (*y(%dS(/3))  -  ^Vds(P)}  ae  +  ... 

{iy(xd  +Vd  +  zmaxs2W)  ~  Vd]  =  0 
-  \c2(P)^j  |  °e  +  {i2z(2yds((3)  +  xdc(p))}  ae  +  ... 

z(Xd  +  Vd  +  ZmaxS  (V7))  —  =  0  (22)-(24) 

These  equations  can  be  solved  by  way  of  the  quadratic  equation.  Examples  show  that  they  will  share  at  least  one 
common,  real  solution  unless  the  coefficients  go  to  zero  (e.g.  xd  =  0  and  (3  =  0  or  180°,  or  when  yd  =  xd  =  0).  When  ix 
=  0,  then  a e  =  2xd/c((3)  .  If  iy  =  0,  then  a e  =  -yd/s((3) .  If  iz  =  0  than  no  useful  information  results. 

A  more  intuitive  understanding  of  these  equations  may  result  when  one  assumes  that  the  deputy  and  deputy  or  in 
the  same  orbit  plane. (i.e.  zmax  =  0).  For  example,  given  the  following  ROE’s  and  associated  LOS  vectors: 

ae  =  unknown  ( but  actually  .04  km) 

xd  =  .01  km 

Vd  =  -.01  km  ix  =  -0.064502287668193 

/ 3  =  557r/ 180  radians  iy  =  0.997917559162865 


Eqns  (22)  and  (23)  yield: 


—0.0791135a;;  +  0.0056437ae  -  .000099  =  0 
0.0791 135a;;  -  0.0056437ae  +  .000099  =  0 

Equation  (24)  has  only  zeros  as  coefficients.  The  allowable  solutions  for  a e  are  0.0313372  or  0.04  km  as  shown  in 
Figure  7. 
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Local  hori/onal  (yaxis)  (km) 


Figure  7.  Possible  solutions  for  ae  given  a  2D  LOS  vector:  ae  =  0.0313372  or  0.04  km. 


The  allowable  solutions  for  ae  result  from  the  fact  that  Eqns  (22)  and  (23)  are  insensitive  to  the  sign  of  the  x  and  y 
components  of  the  LOS  vector  in  the  2D  case.  In  the  case  where  zmax  is  non-zero,  the  solutions  for  Eqns  (22)-(24) 
may  also  result  in  multiple  solutions,  but  they  will  be  skewed  from  the  results  seen  above  based  on  the  values  of  zmax 
and  \| /.  They  will  also  share  only  one  common  solution,  making  the  problem  explicitly  solvable. 

In  the  case  where  xd  is  the  only  unknown,  the  following  quadratic  equations  may  be  derived  from  Eqn  (21): 


{**  -  !}  xd  +  K1  -  il)aec(P)}xd  . . . 

3 


+  {i 


ae(!  -  4C  (P))  +  2 aeyds{(3)  +  yd  +  zmaxs  {$) 


le  2 


c\0)  =  o 


iV(l  -  -c2m  +  (i?  -  1)(2 aeyds((3)  +  y2)  +  i2z2maxs2 &)  -  a2es2(P )  \  =  0 


{i2y  }  x2d  +  {-*Jaec(/3)}  xd  . . . 

3 

\^JJ  1  \"y 

{i2z}  x2d  +  {-i2zaec((3)}  xd  . . . 

•2  I  3 


\  (  (1  A  C  (/^))ae  +  Ud  +  ZrnaxS  (VO  )  zmaxS(^)  }  ~  0 


(25)-(27) 


Note  that  if  ix=0  then  xd  =  ae/2c(P).  If  iy  or  iz  =  0  than  no  useful  information  about  xd  can  be  derived.  Unlike  ae,  all 
three  equations  will  yield  the  same  quadratic  even  in  the  3D  case.  An  example  solution  to  these  equations,  given  the 
ROE's  and  LOS  measurements  below,  is  xd  =  0.00454428  or  0.01  km  as  shown  in  Figure  8. 
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ae  =  .03  km 

Xd  =  unknown  ( but  actually  .01  km) 
yd  =  —.01  km 
(3  =  6l7r/180  radians 


Zmax  —  0.01  km 

•0  =  207t/  180  radians 
ix  =  0.162202680419542 
iy  =  0.965572711469072 
iz  =  0.203370669790363 


Unlike  the  case  with  ae,  the  solutions  to  these  quadratic  equations  are  insensitive  to  the  sign  of  the  x  and  y 
components  of  the  LOS  vector  even  in  the  3D  case.  Other  methods  must  be  applied  to  refine  the  answer.  The 
simplest  method  is  to  see  which  solution  matches  the  sign  on  ix.  To  solve  explicitly,  another  measurement  must  be 
taken.  This  is  explored  in  section  E. 
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Solving  Eq  (21)  for  yd  yields  the  following  quadratic  equations: 


{i-Dvd  +  (2 i2xaes((3)}  yd . . . 

+  | (**  -  l)«d  +  (l  -  il)aec(P)xd  +  ilal(l  -  ^c2(/3))  +  ilz^s2^)  -  ^c2(/3)  | 

{il~1}yd  +  {(*y  -  1)2 aes((3)}  yd  . . . 

+  | ifad  ~  aec((3))xd  +  i2ya2e(  1  -  ^c2(/?))  +  *J4iax«2(V0  -  a2es2(p)\  =  0 
{i2z}yd  +  {i2z2aes{(3)}  yd  . . . 

+  |*2  (ae(!  -  \°2(P))  ~  aeXdC(P)  +X2d  +  zl^82(ll>)j  ~  =  0 

(28)-(30) 

Like  xd,  these  three  equations  will  yield  the  same  results  even  in  the  3D  case.  An  example  solution  to  these 
equations,  given  the  ROE’s  and  LOS  vector  below,  is  yd  =-0.03697245  or  0.03  km  as  shown  in  Figure  9. 

Zmax  =  0.01  km 

ae  —  .04  km  ?/>  =  20-7r/180  radians 

xd  =  .01  km  ix  =  -0.282789522716489 

yd  =  unknown  ( but  actually  .03  km)  %y  =  0.954217664321781 
(3  =  5tt/180  radians  %z  =  0.097461453601231 


Figure  9.  Possible  solutions  for  yd  given  a  3D  LOS  vector:  yd  =-0.03697245  or  0.03  km 
Similarly  to  xd,  the  solution  to  these  quadratic  equations  are  insensitive  to  the  sign  of  the  x  and  y  components  of  the 
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LOS  vector  even  in  the  3D  case.  Other  methods  must  be  applied  to  refine  the  answer.  The  simplest  method  is  to  see 
which  solution  matches  the  sign  on  iy.  To  solve  explicitly,  another  measurement  must  be  taken.  This  is  done  in 
section  E. 


Solving  Eq  (21)  for  zmax  results  in  the  following  equations: 

{*x«2W}  Z?nax  +  ... 

{il  (a2e{^  -  |C2(/3))  +  ae(2yds((3)  -  xdc(0))  +  x2d  +  -  ^c2((5)  +  aexec{(3 )  -  x 

{i2yS2(ip)}  Zliax  +  ... 

Ul  (aei1  ~  \c2m  +  ae(2yds(f3)  -  xdc{(5))  +x2d  +  yj^j  -  a2es2{(5 )  -  2  aeyds((3)  - 
{{i2z  -  1  )s{^)}z2max  +  ... 

| i2z  (a2e{l  -  ^c2(/5))  +  ae{2yds{(3)  -  xdc((3 ))  +xd  +  yd^j  j  =  0 

(31)-(33) 

±V-4  AC 

Solving  using  the  quadratic  equation  results  in  solutions  in  the  form:  2  A  .  Thus  it  can  be  seen  that 

analytically,  the  positive  real  result  is  the  correct  solution. 


Solve  for  \| /  using  the  following  equation 


^  = 


cos  X(1  —  2 K) 


where 


K  = 


-i2x  [a2(l  -  f  c2(/5))  +  ae(2yds((3)  -  xdc(/3 ))  +  x\  +  yj]  +  c2((3 )  -  aexec((3 )  +  x 


j2  p/2 

°xmax 


-i2y  [a2e{  1  -  |c2(/?))  +  ae(2yds{(5)  -  xdc(f3 ))  +  +  2/J]  +  a2s2(/3)  +  2 aeyds(f3)  + 

j2~2 

ymax 


or 


or 


_  -i2  [a2(l  -  f  c2(/3))  +  ae(2ydg(/3)  -  a;dc(/3))  +  + 

(*2  - 

(34)-(36) 

These  solutions  will  only  work  if  ix  or  iy  or  iz  ^  0  respectively,  otherwise  there  will  be  a  zero  in  the  denominator  of 
K.  The  solutions  are  also  unable  to  distinguish  between  quadrants  I  and  II,  or  III  and  IV  because  an  identity  using 
sin(\| /)  was  used  to  derive  the  solutions.  Comparing  the  possible  solutions  to  the  actual  LOS  measurements  is  a 
simple  way  to  resolve  this  problem. 


E.  Analytical  Solutions  with  Two  LOS  Measurements 

Given  two  measurements,  xd  can  be  solved  for  without  the  multiple  solution  problem  found  with  a  single 
measurement.  Solving  equations  9  and  10  for  xd  results  in  the  quadratic  equations  below: 
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Using  the  ix  component  of  the  LOS  vector: 


(4  -  1}  xd  +  {(1  -  **o  Kc(A>)}  xd... 


+  ^  4o 


ae(l  —  ^c2(/^o))  +  2aeya!oS(/?o)  +  2/rfO  +  2Lis2(^o) 


a: 


-  fc2(/3o)HO 


•  <  l 


Jx\ 


9 


1  +  -tu2At2 
4 


i  U|  +  {(i  )^e^(/^l )  3^x1^  A^(CE-eS(/?l  ^  “t-  Udo)  }  Xd  •  •  - 


+  4  4l 


ae(l  —  ^c2(/?i))  +  2aeydos((3i )  +  +  zLj;s2(^i) 


^2(A)|=o 


(37X38) 


Using  the  iy  component  of  the  LOS  vector: 

•  {4}  xd  +  {-i2zOaeC(Po)}  Xd  . .  . 


+  \  ho  (  (!  -  4C  (A) )K  +  2ydos(f3o)ae  +  yd0  +  zmaxs  (ip0)  )  -  zmaxs  (V’o)  }  =  0 


4i 


1  +  L2At2 

4 


xd+  {(1  -  i2zl)aec{(3i)  -  3i2zlujAt(aes(Pi)  +  ydo)}xd... 


+  4  4i 


ae(!  -  ^c2(/?i))  +  2aeydos(p1)  +  y2d0  +  z^s2^) 


(39)-(40) 


Using  the  iz  component  of  the  LOS  vector: 
•  Ko}  Xd  +  {~i2yOaeC(Po)}  Xd  .  .  . 

3 

l,yOu'e\±  ~  VUO ))  T  V‘j/0 


+  1  4oae(!  -  A2  (A)))  +  (ilo  -  l)(2aeyd0s(/30 )  +  j4,)  +  *«o4axs2(^o)  -  a2es2((30)  }  =  0 


•  <  i 


ly  i 
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1  +  -u>2At2 


—  ^-LU2At2  \  X2,  .  .  . 


+  {  — i2i  [aec(/?i)  +  3u>At(ydo  +  aes(/3i))]  +  3u)At(ydo  +  aes(Ai))}  xd  . . . 
3 


+  <i 


4 


ae(l  -  7C  (Ai))  +  2aej/dos(Ai)  +  2/do  +  2  s  (V’l) 


2aet/dos(/3i)  -  yd0  -  a2es2((3 1) 
(41X42) 


Each  pair  of  equations  shares  only  the  correct  solution.  This  can  be  empirically  understood  from  the  fact  that  the 
alternate  solution  for  the  first  equation  results  exists  because  the  component  LOS  value  has  been  squared.  Only  the 
correct  solution  will  match  the  propagated  yao  in  the  second  equation. 
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Given  two  measurements,  yd  can  be  solved  for  without  the  multiple  solution  problem  found  with  a  single 
measurement.  Solving  equations  2  and  3  for  yd  results  in  the  quadratic  equations  below: 


Using  the  ix  component  of  the  LOS  vector: 

•  f2o}  ylo  +  {2^0aes(A)  )}yd  •  •  • 

f  3  a 2 

+  |  (** 0  -  !)*d  +  (!  -  ilo)aec{f5o)xd  +  i2x0a2e(l  -  - c2((30 ))  +  -  -fc2{Po )  =  0 

•  bxi}  2/do  +  {ill  [2aes(^i)  -  3wxdA*]}  j/d0  •  •  • 
i  9 

+  {*: 


,•2 

Xl 


u>2xdAt2  -  SaeUXdAtsiPi)  +  xd  -  aexdc(f3i )  +  a2  (^1  -  ^c2(/3i)  )  +  ^axs2(^i) 
ai 


-x2d  +  aec(Pi)xd  -  yc2(A)}  =  0 


Using  the  iy  component  of  the  LOS  vector: 

•  {ilo  ~  !}  Vdo  +  {{i2yo  ~  1)2 aes(p0)}yd  ■  •  • 

3 


(43)-(44) 


+  s  i 


'yO 


xd  -  aec((3 0)xd  +  ae(  1  -  -c  (Ao))  +  zmaxs  Oo) 


a2es2(Po )  }>  =  0 


byi  -  1}  2/io  +  |(*Ji  -  !)  [2 aes(Ai)  -  3u>zd Ai] }  yd0  •  •  • 
'9 


(Wi-1) 


u>2xdA  t2  —  3aeL0xdAts(Pi ) 


+  i 


•2 

2/1 


^  -  aec((3i)xd  +  ag(l  -  -c2(/?i))  +  z 


2  < 
max'- 


~  aes  (Ai)}  =  0 

Using  the  iz  component  of  the  LOS  vector: 

•  {ilo}y2do  +  {i2zo2aes(p0)}  yd0  . . . 


(45)-(46) 


+  j**o  ^ae(!  -  |c  (A)))  -  aexdc(p0)  +  xd  +  zmaxs  (xp0)  )  -  zmaxs  (V>o)  1  =  0 
{ih}  Vdo  +  {iz l  [2aes(Ai)  -  3cjxd At] }  yd0  •  •  • 

+  {4r9-2-2'  ‘ 

-2L«2(^i)}  =  0 


-u>  xdA t  -  2>aeu)xdAts{Pi)  +  xd-  aezdc(Ai)  +  ae  (  1  -  -c  (Ai)  }  +  zmaxs  (■j/’i) 


(47)-(48) 


Each  pair  of  equations  shares  only  the  correct  solution.  This  can  be  empirically  understood  from  the  fact  that  the 
alternate  solution  for  the  first  equation  results  exists  because  the  component  LOS  value  has  been  squared.  Only  the 
correct  solution  will  match  the  propagated  yao  in  the  second  equation. 


V.  Conclusion 

As  the  need  increases  for  accurate  relative  satellite  navigation  systems  (and  ones  that  can  operate  given  a  sparse 
amount  of  measurement  data),  the  methods  described  in  this  paper  could  become  quite  useful.  Given  only  relative 


^i) 
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angle  measurements  (i.e.  line-of-sight  vectors)  between  two  satellites  and  a  model  based  only  on  two-body 
dynamics,  the  unobservability  problem  is  inevitable.  However,  this  method  allows  one  to  gain  some  knowledge  of 
the  geometric  aspects  of  the  satellite’s  relative  trajectory.  Such  an  algorithm  could  be  incorporated  into  on-board 
flight  software,  employed  by  mission  operators  on  the  ground,  or  simply  used  as  a  post-processing  tool  for  mission 
data.  Future  efforts  will  involve  expanding  the  algorithm  into  a  filter  technique  that  can  estimate  relative  orbit 
elements  given  a  sizeable  amount  of  angle  measurement  data. 
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